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We calculate the temperature and layer-dependent electronic structure of a 20-layer EuO(100)-film 
using a combination of first-principles and model calculation based on the ferromagnetic Kondo- 
lattice model. The results suggest the existence of a EuO (100) surface state which can lead to a 
surface insulator-metal transition. 

In the recent past many theoretical and experimental research works have been focussed on the intriguing properties 
of rare earth metals and their compounds. Among others, the extraordinary surface magnetic properties of the 
lanthanides [Q, as e.g. an enhanced Curie temperature of the Gd(OOOl) surface compared to that of bulk Gd 
have provoked numerous research activities. Concerning the interplay between electronic structure and exceptional 
magnetic properties at the Gadolinium surface a Gd(OOOl) surface state Q is believed to play a crucial role and 
its temperature dependent behaviour has been discussed intensely, e. g, 

Rare-earth materials are so-called local-moment systems, i.e. the magnetic moment stems from the partially filled 
.^/-shell of the rare-earth atom being strictly localized at the ion site. Thus the magnetic properties of these materials 
are determined by the localized magnetic moments. On the other hand, the electronic properties like electrical 
conductivity are borne by itinerant electrons in rather broad conduction bands, e.g. 6s, 5d for Gd. Many of the 
characteristics of local-moment systems can be attributed to a correlation between the localized moments and the 
itinerant conduction electrons. For this situation the ferromagnetic Kondo-lattice model (FKLM) which is also referred 
to as the s-f model has been proven to be an adequate description. In this model, the correlation between localized 
moments and conduction electrons is represented by an intraatomic exchange interaction. 

In this work we will introduce a multi-band FKLM and use it to calculate the temperature and layer-dependent 
bandstructure of EuO films. For all the calculations, the surfaces of the films are parallel to the fcc(lOO) crystal plane 
and the films will be referred to as EuO(lOO) films. These films consist of n equivalent parallel layers. The lattice 
sites within the film are indicated by a greek letter a, /3, 7, . . . , denoting the layer, and by a latin letter i, j, k, . . . , 
numbering the sites within a given layer. The different subbands of the conduction bands of EuO will be denoted by 
the indices m, m' . 

The Hamiltonian for the multi-band FKLM, 

n = ns + nf + nsf, (i) 

consists of three parts. The first, 

— '^ijap '^tama'^jPm'tT, (2) 

contains the kinetic energy of the conduction band electrons, cf^^^ {ciama) is the creation (annihilation) operator of 
an electron with spin a from the m-th subband at the lattice site Rio,. The describe the hopping between the 

m-th subband at the lattice site and the m'-th subband at the lattice site Rj/3. These hopping integrals have to 
be determined within a first-principles band-structure calculation. 

The second part of the Hamiltonian represents the system of the localized /-moments and consists itself of two 
parts, 

^/ = - E J^f^^^^^P -DoJ2 (S-S , (3) 

where the first is the well-known Heisenberg interaction. Here the S^q are the spin operators of the localized magnetic 
moments, which are coupled by the exchange integrals, J^f ■ The second contribution is a single-ion anisotropy term 
which arises from the necessity of having a collective magnetic order at finite temperatures, T > The according 
anisotropy constant Do is typically smaller by some orders of magnitude than the Heisenberg exchange integrals, 
Do « Jtf. 

In addition to the contribution of the itinerant-electron system and the contribution of the localized /-moments we 
have a third contribution to account for an intraatomic interaction between the conduction electrons and the localized 
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/-spins. The form of that contribution in the case of multiple conduction bands can be derived from the general form of 
the on-site Coulomb interaction between electrons of different subbands, with the restriction that electron scattering 
processes caused by the Coulomb interaction are restricted to two involved subbands Q. By distinguishing between 
bands of localized and itinerant electrons, respectively, one gets the on-site Coulomb interaction between both groups 
of electrons. For the special situation of europium with an exactly half-filled localized -.^/-shell and by assuming the 
inter-band exchange to be independent on the different subbands, one eventually gets the multi-band s-/ interaction: 

'Hsf = ——'^^(Tiam ■ Sia- (4) 

iam 

where Sia is the local moment originating from the Hund's rule coupling of the localized /-electrons, (Jiam is the 
Pauli spin operator of an electron from the m-th band of itinerant electrons, both at the lattice site Rm, and J is the 
intraatomic s-/ exchange interaction. 

The many-body problem that arises with the Hamiltonian is far from being trivial and a full solution is lacking 
even for the case of the bulk. In a previous paper we have presented an approximate treatment of the special case 
of a single electron in an otherwise empty conduction band Q . The presented approach applies to the situation in a 
model-type ferromagnetic semiconductor film containing a single-nondegenerated s-like band. However, the approach 
as well as the obtained results are directly transferable to the case of the multi-band Kondo lattice model applicable 
to a real ferromagnetic semiconductor film like EuO. 

Due to the empty conduction bands we are considering throughout the whole paper, the Hamiltonian can be 
split into an electronic part, Tis + 7is/, and a magnetic part, Tif , which can be solved separately For the magnetic 
subsystem Ti f we employ an approach based on the random phase approximation (RPA) which has been described in 
detail in . As the result one gets the layer- and temperature dependent magnetizations of the local-moment system. 
These are depicted for EuO(lOO) films of various thicknesses in Fig. |l|. For the respective calculations, the Heisenberg 
interactions have been approximated within the tight-binding approximation, taking into account the nearest and 
the next nearest neighbour interactions, with values for EuO of Ji/ks — 0.625 K and J2/fcs = 0.125 K, respectively 
The calculated Curie temperature converges as a function of film thickness in a manner similar to that found 
experimentally for Gd(OOOl) films [|lO|. Moreover, the calculated Curie temperature of a 20-layer EuO(lOO) film of 
66.7 K approximates very well the experimental value of Tq for bulk EuO of 69.33 K [pl]| . 

For the electronic part a moment-conserving decoupling approximation (MCDA) for suitably defined Green functions 
is employed ||^ . In the set of equations which constitute the solution for the electronic subsystem, there appear /-spin 
correlation functions which have to be evaluated within the magnetic subsystem (Eq. (|^). These /-spin correlation 
functions mediate the temperature dependence of the electronic structure of the system. 

For the limiting case of ferromagnetic saturation of the localized /-spin system (T = 0), there exists an exact solution 
for the FKLM with empty conduction band {n — 0) While the spin-| spectrum exhibits strong correlation 

effects, the spectrum of the spin-t electron is only rigidly shifted compared to the free solution obtained for the case of 
vanishing s-/ interaction, J = 0. This means, that the spin-l spectra obtained within a ferromagnetic band-structure 
calculation provide the single-particle input (cf. Eq. (|^)) for the FKLM model calculation. Since the solution of the 
FKLM for T = and n = is exact, the problem of double counting of relevant interactions, which would normally 
occur when superimposing the results of a band-structure calculation with a model calculation, is elegantly avoided 
in the presented approach. 

To obtain the temperature-dependent band structure of a EuO(lOO) film, the T — band structures are needed 
as single-particle input for the model calculations. These have been calculated for the Eu-5d bands, which in EuO 



dominate the conduction band, using a standard TB-LMTO program |13 . To account for the film geometry, we 
defined a supercell consisting of 20 EuO(lOO) layers followed by 5 layers of empty spheres, all equidistantly spaced. 
The lattice constant for the EuO(lOO) film has been choosen to be the lattice constant of bulk EuO, a = 5.142 A 
In this approach we neglect surface relaxation and reconstruction, which might be present in EuO films. 

Fig. H shows the temperature- and spin-dependent local densities of states of the first, the second and the center 
layer of a 20-layer EuO(lOO) film with an s-/ exchange interaction of J = 0.25 eV. Here, as in Figs. I and I the energy 
zero refers to the Fermi energy of the system [|l^. For T = 0, the spin-f and the spin- J, spectra occupy the lowest and 
the highest energy range, respectively. With increasing temperature, the spectra for the two spin-directions gradually 
approach each other until for T = Tc {{Sz)/S = 0, fat line) the spectra for both spin-directions become the same. As 
a result, the lower band edge of the spin-f bands is shifted towards lower energies when decreasing the temperature 
below Tc- This effect is known as the red shift of the optical absorption edge in bulk EuO ||ll[. The spin-f densities 
of state for r = {{Sz)/S = 1) are, except for constant energy shift, those obtained within the band-structure 
calculations, which is due to the lack of possibility for the spin-t electron to exchange its spin with the system of 
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perfectly aligned local moments. However, the spin-| electron can flip its spin even for T = by emitting a magnon 
and forming a polaron-like quasiparticle Hence, the spin-| spectra for T = already exhibit correlation effects 
and differ from those of the band-structure calculations. 

In Fig. H the lower band edges of the local densities of states of the surface layer (a = 1) lie at lower energies than 
the respective local densities of state of the center layers and are even touching the Fermi energy, already indicating 
the existence of a surface state band at the lower edge of the Eu-5d bands and a possible surface insulator- metal 
transition. In the following, these interesting indicative observations shall be investigated in more detail. 

In Fig. ||thc T = local spectral density of the surface (a — 1) and the center (a = 10) layer of a 20-layer EuO(lOO) 
film can be seen for both spin directions. A surface state is a state which exists in the so-called forbidden region 
where no bulk states occur. For the bulk spectral density we take that of the center layer (a = 10) of the 20-layer 
film (see Fig. |^). The forbidden regions lie below and above the bulk bands but also includes "white regions" right 
in the middle of the bulk bands. In this respect, in Fig. |^, clearly the lowest bands of the spin-| and of the spin- J, 
spectral densities of the surface layer (a = 1) around the F-point and the M-point constitute surface state bands. 
These surface states originate from the LSDA calculation. With the s-f model calculation one can now investigate 
the temperature dependence of the surface state. 

Fig. U shows the spectral density of the surface layer of a 20-layer film at k = F and k = M for different temperatures. 
For comparison, the respective spectral densities of the center layer of the 20-layer films are plotted. The plotted 
spectral densities clearly indicate surface state bands. For T — {{Sz)/S = 1) the lower band edge of these surface 
state bands lies at the F-point about 0.8 eV and at the M-point about 0.45 eV below the lower conduction-band edge 
of the "bulk" bands of the central layer (cf. dashed line in Fig. |4|). These splittings between the surface states and 
the lower edges of the bulk band are almost independent on the temperature. With increasing temperature both the 
surface states and the lower band edges of the bulk bands for the two spin directions converge Stoner-like. 

As a result, the surface state does not change the down shift of the lower band edge upon cooling below 
Tq (red shift) of the system. However, the lower band edge of the LDOS at the surface is lowered by about 
0.8 eV compared to the bulk-like LDOS in the center of the film. Thus, for the 20-layer EuO(lOO) film, the 
band gap between the occupied 4f^ bands and the unoccupied Sdt^^ bands will be reduced by the same 0.8 eV. 
When decreasing the temperature to T = the red shift will reduce the band gap further: by 0.35 eV accord- 
ing to our calculations and by 0.27 eV according to the experimental results The overall reduction of the 
4f-5dt2g gap will amount to 1.15 eV and 1.07 eV, respectively. These values are exactly in the range of the experi- 
mental 4f-5dt2g gap of bulk EuO at 300 K of 1.12 eV ||ll|. These results indicate a possible surface insulator-metal 
transition in EuO(lOO) films as a function of decreasing temperature, T — > (cf. Fig. ^). Due to the exchange splitting 
of the conduction bands into the lower-energetic spin-f band and the higher- energetic spin- J, band the respective phase 
for T would be a half-metal. As a result, the resistivity of the EuO(lOO) films should be highly dependent on an 
applied magnetic field and a colossal magnetoresistance (CMR) effect would be observed. 
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FIG. 1: Layer-dependent magnetizations, {8^}, of EuO(lOO) films as a function of temperature for Ji/ks = 0.625 K, J2/A;b = 
0.125 K, Do/ks = 0.05 K, for various thicknesses n of the films. Inset: Curie temperatures as a function of film thickness. 
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FIG. 2: Local densities of states of the empty Eu-5ii bands of the first (a = 1), second (a = 2) and center (q — 10) layer of a 
20-layer EuO(lOO) film for different temperatures (Tc = 66.7 K). The energy zero refers to the Fermi energy [|14|. 
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FIG. 3: Local spectral densities S^^ of the surface (a = 1) and the center layer (a = 10) of a 20-layer EuO(lOO) film for T = 
and J = 0.25 eV. The dashed lines in the spectral densities for the center layers (right hand side) reproduce the positions of 
the lower band edges of the surface state bands as given by the spectral densities of the surface layers (left hand side) . The 
energy zero refers to the Fermi energy. 




FIG. 4: Local spectral densities S^^ of the surface (a = 1, solid lines) and the center layer (a = 10, dotted lines) of a 20-layer 
EuO(lOO) film at the F-point and the M-point for J = 0.25eV and for different temperatures (Tc = 66. 7K). The energy zero 
refers to the Fermi energy. 




FIG. 5: Schematics of the surface (half )metal-insulator-transition. For T > Tc the Fermi energy (£f) lies between the 4/^-levels 
and the surface state. For T < Tc E-p lies within the spin-| surface state. 



